Emergence of quasi-condensates of hard-core bosons at finite momentum 
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An exact treatment of the non- equilibrium dynamics of hard-core bosons on one-dimensional lat- 
tices shows that, starting from a pure Fock state, quasi-long-range correlations develop dynamically, 
and that they lead to the formation of quasi-condensates at finite momenta. Scaling relations char- 
acterizing the quasi-condensate and the dynamics of its formation are obtained. The relevance of 
our findings for atom lasers with full control of the wave-length by means of a lattice is discussed. 
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The study of the non-equilibrium dynamics of quan- 
tum gases has been shown to be a very useful tool for the 
understanding of their properties [ij. Since the achieve- 
ment of Bose-Einstein condensation (BEC) 0, where the 
anisotropic expansion of the gas revealed the importance 
of the inter-particle interaction, many experiments study- 
ing the dynamics of such systems have been developed. 
Recently, the possibility of realizing one-dimensional (1- 
D) quantum gases has attracted a lot of attention. They 
can be obtained experimentally in very elongated traps 
[3, or using optical lattices fl. In these quasi- ID systems 
it is possible to go from the weakly interacting regime 
to an impenetrable gas of bosons, i.e. hard-core bosons 
(HCB) (sj. Theoretical studies have analyzed the density 
evolution between these two extreme regimes , showing 
that in general it does not follow a self-similar solution. 
In addition, in the HCB limit the Fermi-Bose mapping 
was generalized to the time-dependent case, where 
the density dynamics revealed dark soliton structures, 
breakdown of the time-dependent mean-field theory, and 
interference patterns of the thermal gas on a ring [Sj. 



In ID systems in equilibrium it was shown that quasi- 
condensates of HCB develop in the homogeneous 0| , har- 
monically trapp ed llOl , and in general cases in the pres- 
ence of a lattice . This is because the ground state of 
these systems exhibits off-diagonal quasi-long-range or- 
der determined by a power-law decay of the one-particle 
density matrix (OPDM) py \xi — la. Illl . i.e. 

there is no BEC in the thermodynamic limit |l2l |. The 
occupation of the lowest natural orbital (NO) (the high- 
est occupied one) is then proportional to y/Nl {Nb being 
the number of HCB) HESIIIl- The NO (<^^) are ef- 
fective single particle states defined as the eigenstates of 
the OPDM [l3j, i.e. they are obtained by the equation 
J2j Pij4>j = : with A,, being their occupations. 

In the present work we study within an exact approach 
the non-equilibrium dynamics of HCB in ID configura- 
tions with an underlying lattice. The presence of the 
lattice enhances correlations between particles, and its 
experimental realization (the so-called optical lattice) al- 
lowed the study of the Mott superfluid-insulator transi- 
tion for soft-core bosons llJI. In contrast to the contin- 



uous case, in a lattice it is possible to create pure Fock 
states of HCB (a HCB per lattice site) where there is 
no coherence in the system. We show that quasi-long- 
range correlations develop in the equal-time-one-particle 
density matrix (ETOPDM) when such states are allowed 
to evolve freely, and that they lead to the formation of 
quasi-condensates of HCB at finite momentum. In ad- 
dition we obtain an universal power law describing the 
population of the quasi-condensate as a function of time, 
independent of the initial number of particles in the Fock- 
state. Finally, we analyze how such systems can be used 
to create atom lasers with a wave-length that can be 
controlled through the lattice parameter. Our analysis is 
based on a generalization of the method used to study 
the ground-state properties of ID HCB on a lattice 
The HCB Hamiltonian on a lattice can be written as 



rii, 



(1) 



with the addition of the on-site constraints for the cre- 
ation (6|) and annihilation (6^) operators: = 6^ = 0, 

and = 1. The hopping parameter is denoted by 

t and Hi = blb^ is the particle number operator. The last 
term in Eq. describes an arbitrary confining potential. 
The Jordan- Wigner transformation , 



i-l i-l 
13=1 13=1 

is used to map the HCB Hamiltonian into the one of 
non-interacting fermions 



XX ni , (3) 



where fj and /j are the creation and annihilation opera- 
tors for spinless fermions, and n{ = //./;• 

This mapping allows to express the equal-time Green's 
function for the HCB in a non-equilibrium system as 

Gy(T) = {^HCB{r)\hb]\^HCB{r)) 

13=1 7=1 



(4) 
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where T is the real time variable. I^'^p^ (''')) the time 
evolving wave- function for the HCB and |^'^(t)) is the 
corresponding one for the non-interacting fermions. 

In what follows, we study the time evolution of ini- 
tial Fock states of HCB once they are allowed to evolve 
freely. Experimentally such states can be created by a 
strong confining potential Va{Ni,a/2)°' ^ t (a is the lat- 
tice constant), which avoid vacancies. Since the equiva- 
lent fermionic system is a non-interacting one, the time 
evolution of an initial wave-fimction (|^'|-)) can be easily 
calculated 



Nj N 



l*F(r)) 



(5=1 <T=1 



which is a product of single particle states, with Nf being 
the number of fermions {Nj = Nb), N the number of lat- 
tice sites and P(t) the matrix of components of 
From here on the method is identical to the one used 
in Ref. the action of 07=1 e~™-^"'-^'' on the right in 
Eq. Q generates only a change of sign on the elements 
PasiT) for cr < .7 — 1, and the further creation of a parti- 
cle at site j implies the addition of one column to P(t) 
with the element PjNf+i{T) = 1 and all the others equal 



to zero (the same applies to the action of 11/3=1 



on the left of Eq. JU). Then the HCB Green's function 
can be calculated exactly as 

Nf + l N Nf + l N 

G^jir) = (0| n E ^^M//' n E ^7^(^)4 io) 

5=1 p=l a=l 7=1 

;p'^(r))'p'«(r)l, (6) 



= det 



where P' (r) and P' (r) are the new matrices obtained 
from P(t) when the required signs are changed and 
the new columns added. The evaluation of Gij (r) is 
done numerically, and the ETOPDM (pijir)) is deter- 
mined by the expression Pijir) ~ (i'lbj'^ ~ Gji{T) + 

5ij{l-2G^^{T)). 

Figure n shows density profiles (a), and their corre- 
sponding momentum distribution functions (MDF) (b) 
for the time evolution of an initial Fock state. Initially, 
the MDF is flat as corresponds to a pure Fock state, and 
during the evolution of the system sharp peaks appear 
at fc = ±7r/2a. Notice that although in the equivalent 
fermionic system the density profiles are equal to the ones 
of the HCB, the MDF remains flat since the fermions do 
not interact, making evident the non-trivial differences 
in the off-diagonal correlations between both systems. 

Since the peaks in the MDF may correspond to 
quasi-condensates at finite momenta, we diagonalize the 
ETOPDM to study the NO. In Fig. E^a) we show the 
lowest NO since is complex) corresponding to 

the results in Fig. ^ for t > (at r = the NO are 
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FIG. 1: (color online). Evolution of density (a) and momen- 
tum (b) profiles of 300 HCB in 1000 lattice sites. The times 
are r = (A), 50ft/t (Q), lOOS/t ( ■ ), and 150fi/i (V). Po- 
sitions (a) and momenta (b) are normalized by the lattice 
constant a. 



delta functions at the occupied sites). They are two- 
fold degenerate, corresponding to right- (for a; > 0) and 
left- (for a; < 0) moving solutions. It can be seen that 
some time after the n = 1 plateau disappears from the 
system, the NO lobes almost do not change their form 
and size (see inset in Fig. Efa) for the right lobe at three 
different times). Furthermore, they move with constant 
velocities vmo = i2at/h (see inset in Fig. |2Ib) for the 
positive one). These are in fact, the group velocities 
V = 1/h dck/dk for a dispersion = —2t cos ka (the 
one of HCB on a lattice) at fc = ±7r/2a. This is further 
confirmed by considering the Fourier transforms of the 
lowest NO Fig. EJb) , which show sharp peaks centered 
at quasi-momenta k = ±7r/2a. The peak at k = +Tr/2a 
appears due to the Fourier transform of the right lobe in 
Fig. 13 (a) (i.e. w |<^0|gi7rx,/2a j^^. ^. > g), and the one 
at k = — 7r/2a due to the fourier transform of the left lobe 
(i.e. w |^0|e-''^^'/2'» for x, <0). Studying the Fourier 
transform of the other NO we have seen that they have 
very small or zero weight at k = ±7r/2a so that we can 
conclude that the peaks at rife=±7r/2a are reflecting the 
large occupation of the lowest NO. 
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FIG. 2: (color online). Lowest NO evolution (a) and its fourier 
transform (b). The times are r = 50h/t (0)i WOh/t (x), and 
150h/t (V). The insets show: (a) the superposed right lobe 
of the lowest NO at r = lOOh/t (x), 125h/t (Q), and 150h/t 
(V); (b) the evolution of the lowest NO right lobe position, 
the line has a slope 2. 
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Whether the large occupation of the lowest NO corre- 
sponds to a quasi-condensate, can be determined study- 
ing the ETOPDM. Fig. O^a) shows the results for \pij\ 
(with i taken at the beginning of the lowest NO left lobe's 
and j > z) at the same times of Figs. ^ El It can be seen 
that off-diagonal quasi-long-range order develops in the 
ETOPDM. It is reflected by a power-law decay of the 
form \pij\ = 0.25 \{xi — Xj)/a|~^/^, that remains almost 
unchanged during the evolution of the system. A care- 
ful inspection shows that this power-law behavior is re- 
stricted to the regions where each lobe of the lowest NO 
exists, outside these regions the ETOPDM decays faster. 
This quasi-coherent behavior in two different segments 
of the system is the reason for the degeneracy found in 
the lowest NO. In addition, a detailed analysis of the 
ETOPDM Fourier's transform shows that the peak in 
the MDF at fc = +Tr/2a is originated by components of 
the ETOPDM with Xi,Xj > 0, and the one at A: = -7r/2a 
by the components with Xi, Xj < 0, so that in the regions 
of the lobes p,j « 0.25 \{x., - ajjO/aj-i/^e^'^^^'-^^^/^a j^^. 
Xi,Xj > 0, and p,j « 0.25 |(a;j - Xj)/a|-i/2e-"(^'-^^)/2a 
for Xi,Xj < 0. The prcfactor of the power law (0.25) was 
found to be independent of the number of particles in the 
initial Fock state. 
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FIG. 3: (color online), (a) ETOPDM for: T^50h/t (Q), 
lOOh/t (x), and 150h/t (V), the line is 0.25 \{xi ~ Xj)/a\^^^^ 
(see text for details), (b) Scaled lowest NO right lobe's vs 
x/Nta for = 101 (Q), 201 ( • ), and 301 (V). The inset in 
(a) shows the maximum occupation of the lowest NO vs Nt 

1/2 

of the initial Fock state, the line is 0.72 A*'^ . 

After having observed off-diagonal quasi-long-range 
correlations in the ETOPDM, it remains to consider how 
the occupation of the lowest NO behaves in the thermo- 
dynamic limit, in order to see whether it corresponds to 
a quasi-condensate. The only new information we need 
at this point is the scaling relation (if any) of the mod- 
ulus of the lowest NO with the number of particles in 
the initial Fock state. For this, we consider the NO at 
those times, where as shown in the inset of Fig. [SJ they 
almost do not change in time, i.e. after the plateau with 
n = 1 disappears. In Fig. Ol^b) we show that a scaled 
NO i\ipo\ = exists when its size is normal- 

ized by Ni,a. (The results for the left lobe arc identical 



due to inversion symmetry.) Figure |2Ib) also shows that 
the lobe size L is approximately N^a. Then evaluating 
Aq = 4>i° Pijff)^ as the double of the integral over a 
single lobe (L ^ a) one obtains 

Xo ^ I dx j dy r°ix)p{x,y)cj)"iy) (7) 

where we did the change of variables x=XNba, y=YNba, 
(tP=N^^ and we notice that the phase factors be- 

tween the NO and the ETOPDM cancel out. The integral 
over X, F is a constant that wc call A. A confirmation of 
the validity of the previous calculation is shown in the in- 
set of Fig.|3Ia). There we plot the maximum occupation 
of the lowest NO (reached when the lobes of the NO have 
a stable form), as a function of the number of particles 
in the initial Fock state. The \/W\j power-law behav- 
ior is evident and a fit allowed to obtain the constant 
A ~ 0.72. Hence, this and the previous results demon- 
strate that the peaks of the MDF in Fig.^b) correspond 
to quasi-condensates with momenta k = -in: /2a. 

The appearance of such quasi-condensates at fc = 
±7r/2a can be understood on the basis of total energy 
conservation. Given the dispersion relation of HCB on 
a lattice, since the initial Fock state has a fiat MDF, its 
total energy is Et = 0. Would all the particles condense 
into one state, it would be the one with an energy cor- 
responding to Zk = Et/N . Taking into account the dis- 
persion relation for HCB on a lattice (e^ = — 2tcos/ca), 
Ik = corresponds to fc = ±7r/2a. Actually, in the ID 
case there is only quasi-condensation, so that the argu- 
ment above applies only to maximize the occupation of 
a given state. In addition, the minimum in the den- 
sity of states at these quasi-momenta reinforces quasi- 
condensation into a single momentum state. 

The process of formation of the quasi-condensate is 
also characterized by a power law. In Fig. ^a) we plot 
the occupation of the lowest NO as a function of the 
evolution time. The log-log scale shows that the occu- 
pation of the quasi-condensate increases quickly and de- 
cays slowly. A more detailed examination shows that 
the population of the quasi-condensate increases in an 
universal way (1.38-\/rt//i, continuous line in Fig. QJa)) 
independently of the initial number of particles in the 
Fock state. The power law is determined by the uni- 
versal behavior of the off-diagonal part of the ETOPDM 
shown before, and the fact that during the formation of 
the quasi-condensate, the NO increases its size linearly 
with time at a rate given by IwatoI = 2at/h. Such a power 
law is followed almost up to the point where the maximal 
occupation is reached. The time at which it is reached 
depends linearly on the number of particles of the ini- 
tial Fock state Tm = Q.32Ni,h/t as shown in the inset of 
Fig. 01 a) . For a tipical experimental setup of rubidium 
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atoms, in a lattice with a recoil energy of Er = 20kHz 
and a depth of 20£'r, the time can be estimated as 
T„ - 5.7iVb(ms). 




FIG. 4: (color online), (a) Time evolution of the lowest NO 
occupation for Nb = 101 (Q), 201 ( ■ ), and 301 (V), the 
straight line is l.SSy rt/h. The inset shows the time at which 
the maximum occupation of the NO is reached (rmt/h) vs Ni,, 
the line is 0.32iV6. (b) Density profiles (inset) and MDF of 
150 HCB evolving only to the right in 1000 lattice sites for 
T = (A), 50ft/t (O), WOh/t ( • ), and 150h/t (V). 

The results above showing that a quasi-coherent mat- 
ter front forms spontaneously from Fock states of HCB, 
suggest that such an arrangement could be used to cre- 
ate atom lasers with a wave-length that can be fully con- 
trolled given the lattice parameter a. No additional effort 
is needed to separate the quasi-coherent part from the 
rest since the quasi-condensate is traveling at the maxi- 
mum possible velocity on the lattice so that the front part 
of the expanding cloud is the quasi-coherent part. The 
actual realization would imply to restrict the evolution of 
the initial Fock state to one direction, as shown in Fig. 
I^b), where we display the MDF for 150 HCB restricted 
to evolve to the right in 1000 lattice sites at the same 
evolution times of Fig. ^ ft can be seen that the values 
of 'T-fc=7r/2(''') sre almost the same in both situations, al- 
though in Fig. Efb) the initial Fock state has almost half 
of the particles. The same occurs to the occupation of 
the lowest NO that in the latter case is not degenerated 
anymore. 

The previous results suggest how to proceed in order 
to obtain lasers in higher dimensional systems where real 
condensation can occur |l6l | . One can employ Mott insu- 
lator states with one particle per lattice site created by a 
very strong on-site repulsive potential U. The latter is re- 
quired in order to obtain a close realization of a pure Fock 
state, since quantum fluctuations of the particle number 
present in a Mott insulator for any finite U ^3 will be 
strongly suppressed. Then the geometry of the lattice 
should be designed in order to restrict the evolution of 
the Mott insulator to one direction only, and to have a 
low density of states around the mean value of energy 
per particle in the initial state. With these conditions 
the sharp features observed in ID should be reproduced 
by a condensate in higher dimensions. 



In summary, we have studied the non-equilibrium dy- 
namics of a Fock states of HCB when they are allowed to 
evolve freely. We have shown that quasi-long-range cor- 
relations develop dynamically in these systems and that 
they lead to the formation of quasi-condensates of HCB 
at quasi-momenta fc±7r/2a. We have studied the dynam- 
ics of the formation of these quasi-condensates, and their 
occupations were found to scale proportionally to ^/Nb. 
These systems can be used to create atom lasers since 
the quasi-condensate develops quickly and decays slowly. 
The wave-length of such lasers is simply determined by 
the lattice parameter. Finally, we have discussed the 
possibility of creating atom lasers in higher dimensional 
systems where true condensation can occur. 
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